3.2.48 \(\int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx\) [148]

3.2.48.1 Optimal result
3.2.48.2 Mathematica [A] (verified)
3.2.48.3 Rubi [A] (verified)
3.2.48.4 Maple [F]
3.2.48.5 Fricas [F]
3.2.48.6 Sympy [F(-1)]
3.2.48.7 Maxima [F]
3.2.48.8 Giac [F]
3.2.48.9 Mupad [F(-1)]

3.2.48.1 Optimal result

Integrand size = 34, antiderivative size = 558 \[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=\frac {(B e-A f) (a+b x)^{1+m} (c+d x)^{1+n} (e+f x)^{-3-m-n}}{(b e-a f) (d e-c f) (3+m+n)}+\frac {(a f (A d f (2+m)+B (d e (1+n)-c f (3+m+n)))+b (B e (d e+c f (1+m))+A f (c f (2+n)-d e (4+m+n)))) (a+b x)^{1+m} (c+d x)^{1+n} (e+f x)^{-2-m-n}}{(b e-a f)^2 (d e-c f)^2 (2+m+n) (3+m+n)}+\frac {((2+m+n) (a b c d f (B e-A f)+b d e ((a B c f+A (b d e-b c f-a d f)) (3+m+n)-(B e-A f) (b c (1+m)+a d (1+n)))-(b c+a d) f ((a B c f+A (b d e-b c f-a d f)) (3+m+n)-(B e-A f) (b c (1+m)+a d (1+n))))-(b c (1+m)+a d (1+n)) (a f (A d f (2+m)+B (d e (1+n)-c f (3+m+n)))+b (B e (d e+c f (1+m))+A f (c f (2+n)-d e (4+m+n))))) (a+b x)^{1+m} (c+d x)^n \left (\frac {(b e-a f) (c+d x)}{(b c-a d) (e+f x)}\right )^{-n} (e+f x)^{-1-m-n} \operatorname {Hypergeometric2F1}\left (1+m,-n,2+m,-\frac {(d e-c f) (a+b x)}{(b c-a d) (e+f x)}\right )}{(b e-a f)^3 (d e-c f)^2 (1+m) (2+m+n) (3+m+n)} \]

output
(-A*f+B*e)*(b*x+a)^(1+m)*(d*x+c)^(1+n)*(f*x+e)^(-3-m-n)/(-a*f+b*e)/(-c*f+d 
*e)/(3+m+n)+(a*f*(A*d*f*(2+m)+B*(d*e*(1+n)-c*f*(3+m+n)))+b*(B*e*(d*e+c*f*( 
1+m))+A*f*(c*f*(2+n)-d*e*(4+m+n))))*(b*x+a)^(1+m)*(d*x+c)^(1+n)*(f*x+e)^(- 
2-m-n)/(-a*f+b*e)^2/(-c*f+d*e)^2/(2+m+n)/(3+m+n)+((2+m+n)*(a*b*c*d*f*(-A*f 
+B*e)+b*d*e*((a*B*c*f+A*(-a*d*f-b*c*f+b*d*e))*(3+m+n)-(-A*f+B*e)*(b*c*(1+m 
)+a*d*(1+n)))-(a*d+b*c)*f*((a*B*c*f+A*(-a*d*f-b*c*f+b*d*e))*(3+m+n)-(-A*f+ 
B*e)*(b*c*(1+m)+a*d*(1+n))))-(b*c*(1+m)+a*d*(1+n))*(a*f*(A*d*f*(2+m)+B*(d* 
e*(1+n)-c*f*(3+m+n)))+b*(B*e*(d*e+c*f*(1+m))+A*f*(c*f*(2+n)-d*e*(4+m+n)))) 
)*(b*x+a)^(1+m)*(d*x+c)^n*(f*x+e)^(-1-m-n)*hypergeom([-n, 1+m],[2+m],-(-c* 
f+d*e)*(b*x+a)/(-a*d+b*c)/(f*x+e))/(-a*f+b*e)^3/(-c*f+d*e)^2/(1+m)/(2+m+n) 
/(3+m+n)/(((-a*f+b*e)*(d*x+c)/(-a*d+b*c)/(f*x+e))^n)
 
3.2.48.2 Mathematica [A] (verified)

Time = 1.09 (sec) , antiderivative size = 508, normalized size of antiderivative = 0.91 \[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=-\frac {(a+b x)^{1+m} (c+d x)^n (e+f x)^{-3-m-n} \left (-((B e-A f) (c+d x))-\frac {(a f (A d f (2+m)+B d e (1+n)-B c f (3+m+n))+b (B e (d e+c f (1+m))+A f (c f (2+n)-d e (4+m+n)))) (c+d x) (e+f x)}{(b e-a f) (d e-c f) (2+m+n)}-\frac {((2+m+n) (a b c d f (B e-A f)-b d e (b (B c e (1+m)+A c f (2+n)-A d e (3+m+n))+a (A d f (2+m)+B d e (1+n)-B c f (3+m+n)))+(b c+a d) f (b (B c e (1+m)+A c f (2+n)-A d e (3+m+n))+a (A d f (2+m)+B d e (1+n)-B c f (3+m+n))))-(b c (1+m)+a d (1+n)) (a f (A d f (2+m)+B d e (1+n)-B c f (3+m+n))+b (B e (d e+c f (1+m))+A f (c f (2+n)-d e (4+m+n))))) \left (\frac {(b e-a f) (c+d x)}{(b c-a d) (e+f x)}\right )^{-n} (e+f x)^2 \operatorname {Hypergeometric2F1}\left (1+m,-n,2+m,\frac {(-d e+c f) (a+b x)}{(b c-a d) (e+f x)}\right )}{(b e-a f)^2 (d e-c f) (1+m) (2+m+n)}\right )}{(b e-a f) (d e-c f) (3+m+n)} \]

input
Integrate[(a + b*x)^m*(A + B*x)*(c + d*x)^n*(e + f*x)^(-4 - m - n),x]
 
output
-(((a + b*x)^(1 + m)*(c + d*x)^n*(e + f*x)^(-3 - m - n)*(-((B*e - A*f)*(c 
+ d*x)) - ((a*f*(A*d*f*(2 + m) + B*d*e*(1 + n) - B*c*f*(3 + m + n)) + b*(B 
*e*(d*e + c*f*(1 + m)) + A*f*(c*f*(2 + n) - d*e*(4 + m + n))))*(c + d*x)*( 
e + f*x))/((b*e - a*f)*(d*e - c*f)*(2 + m + n)) - (((2 + m + n)*(a*b*c*d*f 
*(B*e - A*f) - b*d*e*(b*(B*c*e*(1 + m) + A*c*f*(2 + n) - A*d*e*(3 + m + n) 
) + a*(A*d*f*(2 + m) + B*d*e*(1 + n) - B*c*f*(3 + m + n))) + (b*c + a*d)*f 
*(b*(B*c*e*(1 + m) + A*c*f*(2 + n) - A*d*e*(3 + m + n)) + a*(A*d*f*(2 + m) 
 + B*d*e*(1 + n) - B*c*f*(3 + m + n)))) - (b*c*(1 + m) + a*d*(1 + n))*(a*f 
*(A*d*f*(2 + m) + B*d*e*(1 + n) - B*c*f*(3 + m + n)) + b*(B*e*(d*e + c*f*( 
1 + m)) + A*f*(c*f*(2 + n) - d*e*(4 + m + n)))))*(e + f*x)^2*Hypergeometri 
c2F1[1 + m, -n, 2 + m, ((-(d*e) + c*f)*(a + b*x))/((b*c - a*d)*(e + f*x))] 
)/((b*e - a*f)^2*(d*e - c*f)*(1 + m)*(2 + m + n)*(((b*e - a*f)*(c + d*x))/ 
((b*c - a*d)*(e + f*x)))^n)))/((b*e - a*f)*(d*e - c*f)*(3 + m + n)))
 
3.2.48.3 Rubi [A] (verified)

Time = 0.82 (sec) , antiderivative size = 577, normalized size of antiderivative = 1.03, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.118, Rules used = {172, 172, 27, 142}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (A+B x) (a+b x)^m (c+d x)^n (e+f x)^{-m-n-4} \, dx\)

\(\Big \downarrow \) 172

\(\displaystyle \frac {(a+b x)^{m+1} (B e-A f) (c+d x)^{n+1} (e+f x)^{-m-n-3}}{(m+n+3) (b e-a f) (d e-c f)}-\frac {\int (a+b x)^m (c+d x)^n (e+f x)^{-m-n-3} (b (B c e (m+1)+A c f (n+2)-A d e (m+n+3))+a (A d f (m+2)+B d e (n+1)-B c f (m+n+3))-b d (B e-A f) x)dx}{(m+n+3) (b e-a f) (d e-c f)}\)

\(\Big \downarrow \) 172

\(\displaystyle \frac {(a+b x)^{m+1} (B e-A f) (c+d x)^{n+1} (e+f x)^{-m-n-3}}{(m+n+3) (b e-a f) (d e-c f)}-\frac {-\frac {\int ((m+n+2) (a b c d f (B e-A f)-b d e (b (B c e (m+1)+A c f (n+2)-A d e (m+n+3))+a (A d f (m+2)+B d e (n+1)-B c f (m+n+3)))+(b c+a d) f (b (B c e (m+1)+A c f (n+2)-A d e (m+n+3))+a (A d f (m+2)+B d e (n+1)-B c f (m+n+3))))-(b c (m+1)+a d (n+1)) (a f (A d f (m+2)+B d e (n+1)-B c f (m+n+3))+b (B e (d e+c f (m+1))+A f (c f (n+2)-d e (m+n+4))))) (a+b x)^m (c+d x)^n (e+f x)^{-m-n-2}dx}{(m+n+2) (b e-a f) (d e-c f)}-\frac {(a+b x)^{m+1} (c+d x)^{n+1} (e+f x)^{-m-n-2} (a f (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A f (c f (n+2)-d e (m+n+4))+B e (c f (m+1)+d e)))}{(m+n+2) (b e-a f) (d e-c f)}}{(m+n+3) (b e-a f) (d e-c f)}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {(a+b x)^{m+1} (B e-A f) (c+d x)^{n+1} (e+f x)^{-m-n-3}}{(m+n+3) (b e-a f) (d e-c f)}-\frac {-\frac {((m+n+2) (-b d e (a (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A c f (n+2)-A d e (m+n+3)+B c e (m+1)))+f (a d+b c) (a (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A c f (n+2)-A d e (m+n+3)+B c e (m+1)))+a b c d f (B e-A f))-(a d (n+1)+b c (m+1)) (a f (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A f (c f (n+2)-d e (m+n+4))+B e (c f (m+1)+d e)))) \int (a+b x)^m (c+d x)^n (e+f x)^{-m-n-2}dx}{(m+n+2) (b e-a f) (d e-c f)}-\frac {(a+b x)^{m+1} (c+d x)^{n+1} (e+f x)^{-m-n-2} (a f (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A f (c f (n+2)-d e (m+n+4))+B e (c f (m+1)+d e)))}{(m+n+2) (b e-a f) (d e-c f)}}{(m+n+3) (b e-a f) (d e-c f)}\)

\(\Big \downarrow \) 142

\(\displaystyle \frac {(a+b x)^{m+1} (B e-A f) (c+d x)^{n+1} (e+f x)^{-m-n-3}}{(m+n+3) (b e-a f) (d e-c f)}-\frac {-\frac {(a+b x)^{m+1} (c+d x)^n (e+f x)^{-m-n-1} \left (\frac {(c+d x) (b e-a f)}{(e+f x) (b c-a d)}\right )^{-n} ((m+n+2) (-b d e (a (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A c f (n+2)-A d e (m+n+3)+B c e (m+1)))+f (a d+b c) (a (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A c f (n+2)-A d e (m+n+3)+B c e (m+1)))+a b c d f (B e-A f))-(a d (n+1)+b c (m+1)) (a f (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A f (c f (n+2)-d e (m+n+4))+B e (c f (m+1)+d e)))) \operatorname {Hypergeometric2F1}\left (m+1,-n,m+2,-\frac {(d e-c f) (a+b x)}{(b c-a d) (e+f x)}\right )}{(m+1) (m+n+2) (b e-a f)^2 (d e-c f)}-\frac {(a+b x)^{m+1} (c+d x)^{n+1} (e+f x)^{-m-n-2} (a f (A d f (m+2)-B c f (m+n+3)+B d e (n+1))+b (A f (c f (n+2)-d e (m+n+4))+B e (c f (m+1)+d e)))}{(m+n+2) (b e-a f) (d e-c f)}}{(m+n+3) (b e-a f) (d e-c f)}\)

input
Int[(a + b*x)^m*(A + B*x)*(c + d*x)^n*(e + f*x)^(-4 - m - n),x]
 
output
((B*e - A*f)*(a + b*x)^(1 + m)*(c + d*x)^(1 + n)*(e + f*x)^(-3 - m - n))/( 
(b*e - a*f)*(d*e - c*f)*(3 + m + n)) - (-(((a*f*(A*d*f*(2 + m) + B*d*e*(1 
+ n) - B*c*f*(3 + m + n)) + b*(B*e*(d*e + c*f*(1 + m)) + A*f*(c*f*(2 + n) 
- d*e*(4 + m + n))))*(a + b*x)^(1 + m)*(c + d*x)^(1 + n)*(e + f*x)^(-2 - m 
 - n))/((b*e - a*f)*(d*e - c*f)*(2 + m + n))) - (((2 + m + n)*(a*b*c*d*f*( 
B*e - A*f) - b*d*e*(b*(B*c*e*(1 + m) + A*c*f*(2 + n) - A*d*e*(3 + m + n)) 
+ a*(A*d*f*(2 + m) + B*d*e*(1 + n) - B*c*f*(3 + m + n))) + (b*c + a*d)*f*( 
b*(B*c*e*(1 + m) + A*c*f*(2 + n) - A*d*e*(3 + m + n)) + a*(A*d*f*(2 + m) + 
 B*d*e*(1 + n) - B*c*f*(3 + m + n)))) - (b*c*(1 + m) + a*d*(1 + n))*(a*f*( 
A*d*f*(2 + m) + B*d*e*(1 + n) - B*c*f*(3 + m + n)) + b*(B*e*(d*e + c*f*(1 
+ m)) + A*f*(c*f*(2 + n) - d*e*(4 + m + n)))))*(a + b*x)^(1 + m)*(c + d*x) 
^n*(e + f*x)^(-1 - m - n)*Hypergeometric2F1[1 + m, -n, 2 + m, -(((d*e - c* 
f)*(a + b*x))/((b*c - a*d)*(e + f*x)))])/((b*e - a*f)^2*(d*e - c*f)*(1 + m 
)*(2 + m + n)*(((b*e - a*f)*(c + d*x))/((b*c - a*d)*(e + f*x)))^n))/((b*e 
- a*f)*(d*e - c*f)*(3 + m + n))
 

3.2.48.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 142
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_), x_] :> Simp[((a + b*x)^(m + 1)*(c + d*x)^n*((e + f*x)^(p + 1)/((b*e 
 - a*f)*(m + 1)))*Hypergeometric2F1[m + 1, -n, m + 2, (-(d*e - c*f))*((a + 
b*x)/((b*c - a*d)*(e + f*x)))])/((b*e - a*f)*((c + d*x)/((b*c - a*d)*(e + f 
*x))))^n, x] /; FreeQ[{a, b, c, d, e, f, m, n, p}, x] && EqQ[m + n + p + 2, 
 0] &&  !IntegerQ[n]
 

rule 172
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_)*((g_.) + (h_.)*(x_)), x_] :> With[{mnp = Simplify[m + n + p]}, Simp[ 
(b*g - a*h)*(a + b*x)^(m + 1)*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/((m + 1) 
*(b*c - a*d)*(b*e - a*f))), x] + Simp[1/((m + 1)*(b*c - a*d)*(b*e - a*f)) 
 Int[(a + b*x)^(m + 1)*(c + d*x)^n*(e + f*x)^p*Simp[(a*d*f*g - b*(d*e + c*f 
)*g + b*c*e*h)*(m + 1) - (b*g - a*h)*(d*e*(n + 1) + c*f*(p + 1)) - d*f*(b*g 
 - a*h)*(mnp + 3)*x, x], x], x] /; ILtQ[mnp + 2, 0] && (SumSimplerQ[m, 1] | 
| ( !(NeQ[n, -1] && SumSimplerQ[n, 1]) &&  !(NeQ[p, -1] && SumSimplerQ[p, 1 
])))] /; FreeQ[{a, b, c, d, e, f, g, h, n, p}, x] && NeQ[m, -1]
 
3.2.48.4 Maple [F]

\[\int \left (b x +a \right )^{m} \left (B x +A \right ) \left (d x +c \right )^{n} \left (f x +e \right )^{-4-m -n}d x\]

input
int((b*x+a)^m*(B*x+A)*(d*x+c)^n*(f*x+e)^(-4-m-n),x)
 
output
int((b*x+a)^m*(B*x+A)*(d*x+c)^n*(f*x+e)^(-4-m-n),x)
 
3.2.48.5 Fricas [F]

\[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=\int { {\left (B x + A\right )} {\left (b x + a\right )}^{m} {\left (d x + c\right )}^{n} {\left (f x + e\right )}^{-m - n - 4} \,d x } \]

input
integrate((b*x+a)^m*(B*x+A)*(d*x+c)^n*(f*x+e)^(-4-m-n),x, algorithm="frica 
s")
 
output
integral((B*x + A)*(b*x + a)^m*(d*x + c)^n*(f*x + e)^(-m - n - 4), x)
 
3.2.48.6 Sympy [F(-1)]

Timed out. \[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=\text {Timed out} \]

input
integrate((b*x+a)**m*(B*x+A)*(d*x+c)**n*(f*x+e)**(-4-m-n),x)
 
output
Timed out
 
3.2.48.7 Maxima [F]

\[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=\int { {\left (B x + A\right )} {\left (b x + a\right )}^{m} {\left (d x + c\right )}^{n} {\left (f x + e\right )}^{-m - n - 4} \,d x } \]

input
integrate((b*x+a)^m*(B*x+A)*(d*x+c)^n*(f*x+e)^(-4-m-n),x, algorithm="maxim 
a")
 
output
integrate((B*x + A)*(b*x + a)^m*(d*x + c)^n*(f*x + e)^(-m - n - 4), x)
 
3.2.48.8 Giac [F]

\[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=\int { {\left (B x + A\right )} {\left (b x + a\right )}^{m} {\left (d x + c\right )}^{n} {\left (f x + e\right )}^{-m - n - 4} \,d x } \]

input
integrate((b*x+a)^m*(B*x+A)*(d*x+c)^n*(f*x+e)^(-4-m-n),x, algorithm="giac" 
)
 
output
integrate((B*x + A)*(b*x + a)^m*(d*x + c)^n*(f*x + e)^(-m - n - 4), x)
 
3.2.48.9 Mupad [F(-1)]

Timed out. \[ \int (a+b x)^m (A+B x) (c+d x)^n (e+f x)^{-4-m-n} \, dx=\int \frac {\left (A+B\,x\right )\,{\left (a+b\,x\right )}^m\,{\left (c+d\,x\right )}^n}{{\left (e+f\,x\right )}^{m+n+4}} \,d x \]

input
int(((A + B*x)*(a + b*x)^m*(c + d*x)^n)/(e + f*x)^(m + n + 4),x)
 
output
int(((A + B*x)*(a + b*x)^m*(c + d*x)^n)/(e + f*x)^(m + n + 4), x)